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We investigate the combined effect of spin-orbit scattering and electron-electron interactions on 
the probability distribution of (/-factors of metal nanoparticles. Using random matrix theory, we find 
that even a relatively small interaction strength significantly increases <?-factor fluctuations for not- 
too-strong spin-orbit scattering (ratio of spin-orbit rate and single-electron level spacing 1/t so <5 < 1), 
and leads to the possibility to observe ^-factors larger than two. 
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Electronic properties of metal nanoparticles can be 
studied on the single-electron level using "tunneling spec- 
troscopy" , the measurement of the conductance of a 
metal particle coupled to source and drain electrodes 
via tunneling contacts [1]. These measurements have 
revealed important insights into the nature of the elec- 
tronic ground and of individual excited states of normal 
metal, ferromagnetic, and superconducting nanoparticles 
[2, 3, 4, 5, 6]. Although tunneling spectroscopy involves 
processes in which only a single electron is added to or 
removed from the metal particle, the location of conduc- 
tance peaks gives information about many-electron en- 
ergy levels and, hence, about the role of electron-electron 
interactions. 

The combination of electron-electron interactions and 
mesoscopic fluctuations of the density of states of a 
normal-metal nanoparticle can lead to a many-electron 
ground state with nontrivial spin 5 = 1,5* = 3/2, or even 
larger [7, 8]. (Without interactions, the ground state spin 
is always 5 = or 5 = 1/2.) In principle, the spin of a 
many-electron state |fc) can be measured via the deriva- 
tive of the energy Ek versus an applied magnetic field H , 
which is parameterized using the "g- factor" g k , 
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If spin is a good quantum number, the positions of tun- 
neling spectroscopy conductance peaks are unaffected by 
a nontrivial value of 5, since tunneling spectroscopy mea- 
sures differences of (/-factors of many-electron states for 
which the electron number differs by one. Since the 
spin of a nanoparticle changes by 1/2 upon addition or 
removal of an electron, all observed (differences of) g- 
factors are equal to two, irrespective of the individual 
g-factors of the two many-electron levels participating in 
the transition. 

In this letter, we consider (/-factor differences mea- 
sured in tunneling spectroscopy in the presence of both 



electron-electron interactions and spin orbit scattering. 
When spin and orbital degrees of freedom are coupled, 
randomness in the orbital part of wavefunctions is passed 
on to the spin part. Randomizing the electron spin lifts 
the "selection rule" that prohibited the observation of g- 
factors larger than two in the absence of spin-orbit scat- 
tering. It also leads to a decrease of ^-factors, a sup- 
pression of the long-range exchange interaction (which 
is responsible for the high-spin states), and level-to-level 
fluctuations of (/-factors [9, 10]. As we show here, selec- 
tion rules are lifted already for a small spin-orbit scat- 
tering rate l/r so < 5, S being the mean spacing between 
single-electron levels, whereas the decrease of (/-factors 
and the suppression of the exchange interaction become 
effective at a larger spin-orbit rate l/r so > 5 only, leav- 
ing a substantial parameter window where (/-factor dif- 
ferences larger than two can be observed. Such large 
(/-factors are a true many-electron phenomenon since, 
without interactions, all measured (/-factors correspond 
to single-electron levels and are always < 2 [11]. 

Model. Without interactions, the single-electron wave- 
functions and energy levels of a metal nanoparticle are 
described by random matrix theory. With spin orbit scat- 
tering, the appropriate random matrix ensemble interpo- 
lates between the Gaussian Orthogonal Ensemble (GOE) 
and the Gaussian Symplectic Ensemble (GSE) [9], 
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Writing the spin degrees of freedom explicitly, one has 

A 3 

H GOE = S(g>l 2l H so = -—y^iAj^o-j, 

3 = 1 

where I2 is the 2x2 unit matrix in spin space, (Tj is the 
Pauli matrix (j = 1,2,3), 5 is an N x N real symmet- 
ric matrix, Aj is an N x N real antisymmetric matrix 
(J = 1,2,3), and A 2 = it/t so S is the dimensionless spin- 
orbit scattering rate. The elements of the matrices 5, 
A±, A2, and A3 are drawn from independent Gaussian 
distributions with zero mean and with equal variances 
for the off-diagonal elements. The diagonal elements of 
5 have double variance, whereas the diagonal elements 
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of A\, A 2 , and A3 are zero because of the antisymme- 
try constraint. The limit N — > 00 is taken at the end 
of the calculation. Each eigenvalue e M of H is doubly 
degenerate, with wavefunctions Wi and tp^ related by 
time-reversal. 

In normal-metal nanoparticles, the main contribution 
to electron-electron interactions is described by the "con- 
stant exchange interaction model" [12, 13], 



H ex — — JS , 
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where S is the total spin of the particle. The ratio of 
the exchange constant J and the mean spacing between 
single-electron levels 6 corresponds to one of the Fermi 
Liquid constants of the metal. For most normal metals 
one has 0.2 < J/6 < 0.4, in agreement with electron- 
liquid theory [14] , although smaller and larger values oc- 
cur as well (see the discussion at the end of this letter). 
Combining the constant exchange interaction (3) and the 
single-electron Hamiltonian (2), and including the Zee- 
man coupling to a magnetic field H in the z direction 
[11], one has 



H 



E + ^2) - JS 2 ^ B HS Z , (4) 



where the first term comes form the diagonalization of 
the single-electron Hamiltonian (2). 

If the number of electrons N e is even, all many-electron 
states are non-degenerate in the presence of spin-orbit 
scattering. Hence all even-electron states have (/-factors 
equal to zero. With spin-orbit scattering, odd-electron 
states are twofold degenerate (Kramers' degeneracy). 
Since tunneling spectroscopy measures differences of g- 
factors for many-electron states with N e and N e + 1 elec- 
trons, a (/-factor measured using tunneling spectroscopy 
is the (/-factor of an odd-electron state if spin-orbit scat- 
tering is present. We assume that the nanoparticle re- 
laxes to the even-electron ground state between tunneling 
events so that one measures ^-factors of the odd-electron 
ground state and the odd-electron excited states upon in- 
creasing the bias voltage. In the (generic) case that the 
two tunneling contacts between the nanoparticle and the 
source and drain reservoirs have very different conduc- 
tances, the height of a conductance peak is proportional 
to the matrix element [1] 

w k = ]T |(JV e + l,£#t( r) |7V e , )| 2 , (5) 

<T=U 

where the creation operator V>£(r) creates an electron 
with spin a in the grain at the position r of the point 
contact with the smaller conductance, \N e , 0) is the even- 
electron ground state, and \N e + 1, k) is an odd-electron 
excited state. Without spin-orbit scattering, the weights 
Wk are zero for those states | N e + 1, k) for which the spin 



S z differs by more than 1/2 from the spin of the even- 
electron ground state | N e , 0). (This is the "selection rule" 
referred to in the introduction.) 

Limit of weak spin-orbit scattering. We first address 
the limit A <C 1, for which perturbation theory in A 
is possible. In view of the above comments, we need 
to consider (/-factors of odd-electron states only. For an 
odd-electron state |fc) with spin S = 1/2 without spin- 
orbit scattering, spin-orbit scattering only affects the spin 
contribution to the (/-factor to quadratic order in A [15]. 
However, all nanoparticles have low- lying odd-electron 
states with spin S = 3/2 if the exchange interaction is 
present. (There even is a small but nonzero probabil- 
ity that the odd-electron ground state has spin 5 = 3/2 
[7, 8].) Spin-orbit scattering lifts the fourfold degeneracy 
of an 5 = 3/2 odd-electron state and splits this state into 
two doublets. As we show below, spin-orbit scattering de- 
termines the g-factors of these doublets already to zeroth 
order in A, whereas the matrix elements determining the 
corresponding tunneling spectroscopy peak height Wk are 
nonzero with finite probability. 

Labeling the four members of the S — 3/2 quadruplet 
by the z component of the spin, S z = p — 5/2, p — 
1, 2, 3, 4, the matrix elements of H so can be arranged in 
a 4 x 4 matrix Vk of the form 
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with a and d real numbers and b and c complex num- 
bers. The specific form of (6) follows from time-reversal 
symmetry and guarantees that that the eigenvalues of 
Vk are doubly degenerate, in accordance with Kramers' 
theorem. One has Vk = to first order in H so , since 
spin-orbit scattering does not mix states with opposite 
spin belonging to the same energy level. Calculating 14 
to second order in H so , we consider the special case when 
the 5 = 3/2 quadruplet is split by virtual transitions to 
one nearby odd-electron state \l) only. We refer to Ref. 
16 for a general discussion; neglecting other states is jus- 
tified if the energy difference of the virtual transition we 
consider is much smaller than other energy separations. 

Since we are interested in g-factors only, it is sufficient 
to calculate the ratios b/d and c/d, for which we find 
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where /1 and v refer to the two single-electron levels in- 
volved in the virtual transition. The 5 = 3/2 quadruplet 
splits into two doublets with (/-factors gk\ and gki, with 
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FIG. 1: Ensemble averaged g-factors for J/5 = 0, 0.1, 0.2, 
0.3, 0.4, 0.5, and 0.6. 

g\ 2 = 48-3<&, (8) 

Whether g k i or g k2 corresponds to the lower lying doublet 
depends on the relative position and spin of the virtual 
state; if the unperturbed 5 = 3/2 state is the ground 
state, the doublet with lower energy has (/-factor g k i [16]. 
Using the distribution of the matrices Ai, A2, and A3, 
one finds that g k \ has the distribution 

Pi(9ki) = \-=M=, 2<g kl <4. (9) 
2 ^48 - 3g 2 kl 

The second g-factor g k2 takes values in the interval < 
g k2 < 6 and is related to g k \ via Eq. (8). 

To zeroth order in A, the tunneling spectroscopy peak 
heights (5) are nonzero if the even-electron ground state 
I N e , 0) has spin 5=1. This occurs with significant prob- 
ability for J > 0.3(5 [7, 8], so that the anomalous g- factors 
g k \ and g k 2 are visible in a tunneling spectroscopy exper- 
iment with finite probability. Increasing the spin-orbit 
scattering rate A further increases the visibility of peaks 
with the largest (/-factors by mixing even-electron states 
with 5 = and 5=1. Indeed, for J = 0.3(5 the av- 
erage energy difference between the lowest lying 5 = 
and 5 = 1 even-electron states is <~ S — 2 J = 0.4(5, so 
that even moderate spin-orbit scattering (A > 0.58) has 
a matrix element XS between the two states that is com- 
parable with the energy difference. Inclusion of virtual 
excitations to other excited state changes the formulas 
for g k i and g k 2, but not the conclusion that spin-orbit 
scattering affects the g- factors of the 5 = 3/2 states to 
zeroth order in A. 

Arbitrary spin-orbit scattering rate. In order to ad- 
dress the effects of a finite spin-orbit scattering rate, we 
have numerically diagonalized the Hamiltonian (4) for 
< A < 2.8. We first diagonalized Hq and considered 
the interaction iJ ex in the basis of the 92 (76) lowest lying 
many-electron eigenstates of H n for N e odd (even). We 
then diagonalized the remaining many-electron Hamilto- 
nian and calculated the g-factors g k of the M — 8 lowest- 



lying odd-electron states, together with peak height w k 
for transition from the even-electron ground state, see Eq. 
(5). The random matrices in our simulation are taken of 
size 27V = 400 for 2 < A < 2.8 and of size 2N = 200 for 
A < 2. We averaged over 300 realizations of the random 
matrices, corresponding to a mesoscopic average over an 
ensemble of nanoparticles with equal size and spin-orbit 
scattering rate but different disorder configurations. In 
the analysis of the numerical data, we discarded all levels 
for which the peak height w k is below a threshold w tr , 
which we arbitrarily set at 

w t r = 0.1 x max^Wfc. (10) 

The threshold mimics the experimental reality that small 
peaks cannot be distinguished from the noise, and, hence, 
have their g- factors left out in the statistical analysis. 
Further, omitting g-factors for which w k < w tI enforces 
the "selection rules" in the absence of spin-orbit scatter- 
ing. We verified that the precise definition of w tI does 
not affect our conclusions. 

Interaction effects increase the ensemble averaged g- 
factor (g) significantly for J > 0.38 and A < 2, see Fig. 
1. In fact, there is a substantial parameter window for 
which (g) > 2. By itself, such an increase of (g) has lim- 
ited experimental relevance, since there is no independent 
method to measure A. In fact, comparison of the mea- 
sured (g) with theory is used to determine the spin-orbit 
rate l/r so = \ 2 S/ir in nanoparticles [5, 6]. This problem 
does not exist for the full (cumulative) probability distri- 
bution of g- factors (average and fluctuations), which is 
shown in Fig. 2. The values of J and A in Fig. 2 are cho- 
sen such that all distributions have the same average (g) . 
The distributions of Fig. 2, together with our numerical 
results for other values of J and A (not shown), show 
that the exchange interactions substantially enhance the 
(/-factor fluctuations (at the same value of the average) . 
The probability P(g > 2) to find a (/-factor larger than 
two is shown in the inset of Fig. 2, as a function of J and 
A. 

Presently, (/-factor distributions have been measured 
for the noble metals only [5, 6], for which J/8 < 0.1 and 
interaction effects arc negligible [17]. Indeed, the distri- 
butions measured in Refs. 5, 6 are in good agreement 
with the non-interacting theory [9, 10]. Interaction ef- 
fects, nontrivial spin states, and, hence, g-factors larger 
than two should be observable for most other metals. For 
alkali metals, J/8 is in the range 0.2-0.3 [18], as well for 
Ti, Zr, and Mg [19, 20]. Even stronger interaction effects 
are expected for Nb, Rh, and Y nanoparticles, for which 
spin-density calculations set J/8 around 0.4 [21], and for 
Pt, Pd, and V, which have 0.6 < J/8 < 1 [22, 23]. For 
particles in the nm size range, spin-orbit effects are ex- 
pected to be moderate or weak (dimcnsionless spin-orbit 
rate A of order unity or smaller), except for the elements 
with the highest atomic numbers (Au, Pt). (Measured 
spin-orbit rates in Ag and Cu nanoparticles with radius 
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strong interaction effects. 
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FIG. 2: Cumulative ^-factor distribution for A = 0.70, J = 
(solid curve), A = 0.85, J = 0.25 (points), A = 0.90, J = 0.3(5 
(dotted), A = 1.0,d and J = 0.45 (dash-dot), and A = 1.15, 
J = 0.65 (dashed). The values of A are chosen such that 
(g) — 1.58 in all cases. Inset: Probability for a level to have 
a p-factor larger than two, for J/8 = 0.1 — 0.6, versus A. 



~ 4 nm were in the range A ~ 1 [5].) For those elements 
for which spin-orbit effects are too weak to make large 
g-factors visible, spin-orbit scattering can be enhanced 
by doping with a small amount of, e.g., Au atoms [24]. 
Similarly, doping with ferromagnetic atoms may increase 
J and drive the metals towards the Stoner instability at 
J/S=l [25]. 

In conclusion, we have shown that the presence of 
an exchange interaction with strength J > 0.28 leads 
to a significant broadening of the probability distri- 
bution of tunneling spectroscopy g-factors in normal 
metal nanoparticlcs. In particular, g-factors larger than 
two can be observed, which are a signature of non- 
trivial many-electron states. Weak spin-orbit scattering 
(l/ T so 8) is crucial in rendering the large g- factors 
observable, since it mixes many-electron states with dif- 
ferent spin and, hence, lifts spin selections rules. It is 
only at larger spin-orbit scattering rates l/r so 3> 8 that 
spin-orbit scattering fully randomizes the spin and sup- 
presses Zeeman contribution to g-factors and the matrix 
elements of the exchange interaction. The interaction 
range J/8 > 0.2 is appropriate for most metals, except 
for Al and the noble metals, which have negligible ex- 
change interaction effects. 

Whereas in gated semiconductor quantum dots the 
existence of nontrivial spin states could be inferred in- 
directly from the statistical distribution of Coulomb- 
blockade peak spacings [26], a comparison of the de- 
tailed parameter-dependence of successive peak positions 
[27], or from the succession pattern of (/-factors of many 
consecutive levels [28], no such methods are available in 
metal nanoparticlcs without a gate or with a very lim- 
ited gate voltage range. We hope that our finding that 
spin-orbit scattering renders nontrivial g-factors visible 
in standard tunneling spectroscopy without a gate mo- 
tivates further experiments on metal nanoparticles with 
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